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CHAPTER  I 

INTRODUCTION 

1 . 1  Fundamental  Characteristics  of  Cavitat Ion 

When  the  local  pressure  falls  below  a  critical  value,  small  bubbles 
grow  by  vaporization  or  by  degassing.  Bubble  growth  by  vaporization  is 
called  vaporous  cavitation  if  caused  by  dynamic  pressure,  reduction  at 
essentially  a  constant  temperature.  For  vaporous  cavitation,  which  is 
of  primary  interest,  the  critical  value  of  pressure  is  less  than  or 
equal  to  the  vapor  pressure  of  the  liquid.  In  the  initial  steps  of 
the  phenomenon,  mic.rohubbles  grow  in  the  minimum  pressure  region  and 
then  when  transferred  to  higher  pressure  regions  by  the  fluid  flow, 
the  growth  will  be  arrested.  The  collapse  of  bubbles  will  cause 
material  damage,  noise  and  vibration. 

On  the  other  hand,  a  further  reduction  in  the  minimum  pressure 
will  promote  the  growth  of  bubbles  and  may  result  in  large  cavities 
attached  to  the  body.  The  succeeding  stages  are  distinguished  from 
the  initial  stage  by  the  term  developed  cavitation,  which  may  produce 
a  serious  performance  loss  when  it  occurs  on  hydromachinery,  marine 
vehicles,  etc. 

1 . 2  The  Effect  of  the  Surface  Irregularities  on  Cavitation 
Basically,  surface  Irregularities  may  be  divided  into  two  categories 

namely,  isolated  and  distributed. 

Isolated  irregularities  will  cause  high  local  velocities,  low 


pressures  and  turbulence  in  the  neighborhood  of  the  projection. 


Distributed  irregularities  will  change  velocities,  pressure  and 
urbulence  generally  throughout  the  volume  of  boundary  layer  fluid. 
Surface  irregularities,  therefore,  advance  the  onset  of  cavitation  on 
a  body  and  reduce  its  design  performance.  In  general,  isolated 
irregularities  are  more  harmful  to  cavitation  performance  than 
distributed  irregularities  because  the  local  pressure  reduction  is 
more  severe. 

1. 3  Thesis  Objectives 

The  objective  of  this  thesis  is  to  investigate  the  effect  of 
pressure  gradient  on  the  cavitation  inception  of  an  isolated  irregularity. 
In  so  doing,  a  curved  wall  element  was  attached  on  one  wall  in  the 
water  tunnel  test  section,  and  triangular  protrusions  of  various  sizes 
were  mounted  on  the  opposite  side  as  isolated  irregularities.  In 
order  to  attain  the  objective,  the  velocity  profiles  of  the  boundary 
layer  and  the  pressure  gradient  were  measured  in  the  initial  phase  of 
the  test  and  then  subsequently  the  cavitation  number  was  measured  for 
the  cases  in  which  the  tip  of  the  triangular  protrusion  was  both 
inside  and  outside  the  boundary  layer. 
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CHAPTER  II 

PREVIOUS  INVESTIGATIONS 

2.1  The  Influence  of  Irregularities  on  Cavitation 

Surface  irregularities  are  one  of  the  major  factors  influencing 
cavitation  inception.  Recent  requirements  of  compact  and  high  speed 
hydromachinery  have  demanded  more  extensive  knowledge  of  the  effect 
of  irregularities  on  cavitation. 

The  initial  work  in  this  area  was  concerned  with  performance 
losses  in  turbomachinery.  Shalnev  [1]  studied  cavitation  caused  by 
surface  irregularities  in  the  clearance  between  the  pump  blades  and 
the  casing  wall.  However,  his  results  would  only  be  applicable  to 
similar  boundaries  since  they  were  not  correlated  with  boundary 
layer  parameters. 

Walker  [2]  employed  the  concept  of  relative  height  of  roughness 
(h/6)  ,  where  h  is  the  height  of  roughness  .and  6  the  local  boundary 
layer  thickness.  However,  his  investigation  was  conducted  for  specific 
conditions  and  was  not  done  systematically. 

Holl  [3,4]  systematically  investigated  cavitation  inception 
numbers  for  two  dimensional  isolated  roughness  elements  ,  namely, 
circular  area  and  sharp  edged  triangular  strips  on  a  flat  plate. 

By  placing  the  roughness  at  different  locations  in  the  boundary  layer 
of  a  smooth  wall,  he  found  that  the  cavitation  number  could  be  correlated 
with  the  relative  height  of  roughness  (h/o),  a  Reynolds  number  (Re^) 
based  on  roughness  height  and  local  velocity,  and  the  ordinary  Karman- 
Pohlhausen  type  shape  factor  (H) . 
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Benson  [5]  investigated  cavitation  inception  numbers  for  three- 
dimensional  isolated  roughness  elements,  namely,  hemispheres,  cones. 
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and  circular  cylinders  on  a  flat  plate.  He  showed  that  the  cavitation 
correlated  well  with  a  local  Reynolds  number  based  on  the  roughness 
height  and  the  flow  velocity  at  the  height  of  the  roughness. 

The  effect  of  distributed  roughness  on  cavitation  is  more 
complicated  than  the  case  of  an  isolated  roughness  because  there  is  an 
interrelation  between  the  boundary  layer  development  and  the  degree 
of  surface  roughness.  However,  several  studies  have  been  conducted 
since  the  first  publication  to  appear  regarding  distributed  roughness 
by  Benson  [6] . 

The  classical  work  in  distributed  roughness  is  that  of  Arndt  [7] 
and  Arndt  and  Ippen  [8],  From  an  examination  of  the  results  compiled 
by  Benson  [6},  Arndt  indicated  that  the  cavitation  occurring  in  the 
roughness  grooves  was  probably  of  the  gaseous  type,  and  the  cavitation 
occurring  in  the  boundary  layer  was  vaporous  cavitation.  Arndt  and 
Ippen  studied  the  effect  of  distributed  roughness  by  experiments 
with  triangular  groove  irregularities.  They  showed  that  cavitation 
occurred  in  the  middle  of  the  boundary  layer,  and  correlated  the 
incipient  cavitation  number  (a^)  with  the  wall  shear  stress  coefficient 
( C^)  for  the  rough  boundary. 

Additional  studies  concerned  with  distributed  roughness  on  cavitation 
are  (1)  experiments  with  sand  grain  irregularities  by  Messenger  [9] 
and  Huber  [10];  (2)  experiments  with  wire  screen  irregularities  by 
Bechtel  [11],  and  (3)  the  effect  of  surface  roughness  on  the  performance 
of  hydrofoils  by  Numachi  et  al„  [12,13].  The  investigations  by  Messenger, 
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Huber  and  Bechtel  confirmed  the  linear  relationship  between  and 
which  was  established  by  Arndt  [7]  and  Arndt  and  Ippen  [8], 

2.2  Velocity  Profiles  of  Boundary  layers 

The  velocity  profiles  for  turbulent  boundary  layers  may  be 
divided  into  three  parts: 

1)  Inner  layer 

2)  Outer  layer 

3)  Overlap  layer. 

For  the  inner  layer,  dominated  by  viscous  shear,  Prandtl  [14] 
deduced  that  the  mean  velocity  (u)  must  depend  upon  the  wall  shear 
stress,  the  fluid  properties,  and  the  distance  (y)  from  the  wall. 

Thus,  he  proposed  the  inner  law  or  so-called  law  of  the  wall: 

u  =  f(T  ,  P,  M.  y) ,  (2-1) 

w 

where  T  is  the  shear  stress  at  the  wall, 
w 

For  the  outer  layer,  dominated  by  turbulent  shear.  Von  Karman  [15] 
deduced  that  the  wall  tends  to  act  merely  as  a  source  of  retardation, 
reducing  the  local  velocity  (u)  below  the  Ireestream  value  (U)  in  a 
manner  which  is  independent  of  viscosity  (u)  but  dependent  upon  the 
shear  stress  at  the  wall  and  the  distance  (y)  over  which  its  effect 
has  diffused.  Thus,  he  proposed  the  outer,  or  velocity  defect,  law: 

U  =  u  =  f(T  ,  p,  y,  5) ,  (2-2) 

w 


where  6  is  the  boundary  layer  thickness. 
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For  the  overlap  layer,  Millikan  fl6]  simply  specified  that  the 
two  functions  merge  together  smoothly  over  some  finite  region  between 
the  inner  and  outer  layer. 

From  dimensional  analysis,  the  non-dimensional  form  of  the  law 
of  the  wall  and  the  velocity  defect  law,  respectively,  are 


u 


P’J  y 
f(— > 


and 


where 


(2-3) 


(2-4) 


(2-5) 


In  the  overlap  layer,  a  velocity  profile  can  be  expressed  by 
using  a  mixing  length  model.  The  Reynolds  stress  in  turbulent 
flow  is  expressed  as. 


T 

t 


C  p  i2 


9u  3u 
3y  3y  * 


(2-6) 


where  £  is  a  mixing  length  and  C  is  a  coefficient  of  order  one.  This 
is  called  Prandtl' s  mixing-length  hypothesis.  The  eddies  involved  in 
momentum  transfer  have  characteristic  vorticities  of  order  v' / Z  and  they 
maintain  their  vorticity  because  of  their  interaction  with  the  mean 
shear  (t~)  ,  so  that 
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where  v'  is  the  instantaneous  velocity  in  the  y  direction. 

From  equations  (2-6)  and  (2-7) 

tt  -  p  V’  <l|a  .  (2-8) 

On  the  assumption  that  u'  and  v'  are  well  correlated „  the  Reynold's 
stress  will  be 

Tt  =  pu*\  (2-9) 


where  u'  is  the  instantaneous  velocity  in  the  x  direction  and  uy.  is 
the  wall-friction  velocity. 

From  equations  (2-8)  and  (2-9) 


J_  3u  Jl 
u*  3y  =  k  y  ’ 

where  k  is  a  const ant. 

After  integrating  equation  (2-10), 

u  1  . 

-  —  =  7-  icn  y  +  const. 
u*  k 


(2-10) 


(2-11) 


Furthermore,  equation  (2-11)  can  be  written  as 


u 


where  k  is  known  as  the  constant  of  Von  Karman  and  B  is  an  empirical 
constant.  Nikuradse  [17]  in  1930  found  that  the  constants  k  and  B 
were  0.40  and  5.5  and  Coles  [18]  in  1955  employed  0.41  and  5.0  for 
k  and  B,  respectively. 
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2 . 3  The  Law  of  the  Wake  for  Turbulent  Boundary  Layers 

The  outer  law  is  more  sensitive  to  external  parameters,  particularly 
the  pressure  gradient,  as  shown  in  Fig.  6-4,  p.  470  of  Ref.  [19], 

Clauser  [20,21]  was  the  first  investigator  to  introduce  new 
parameters  involving  the  pressure  gradient  for  the  outer  layer; 
these  are  given  by 


.(.X  A  iE) 

*•^6  *  T  dxj  ’ 
w 


(2-13) 


After  further  investigation,  he  concluded  that  6  should  be  replaced 
by  the  displacement  thickness.  Therfore,  the  accepted  parameter  (3)  is 

'41-  <2-i4> 

w 


A  shape  factor  (G)  which  remains  constant  in  an  equilibrium 
boundary  layer  was  defined  by  Clauser  as 


G 


dy, 


(2-15) 


where 


(2-16) 


Equation  (2-15)  can  be  related  to  the  ordinary  Karman-Pohlhausen  type 
shape  factor  (H)  by 
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The  shape  factor  (H)  is  not  constant  in  an  equilibrium  boundary  layer 
because  the  skin  friction  varies  with  x.  It  is  said  that  G  is  a  better 
shape  factor  than  H  in  terms  of  predicting  the  properties  of  equilibrium 
boundary  layers.  Nash  [22]  approximated  G  by  using  the  pressure 
gradient 

G  ~  6.1  /ITT  OT  -  1.7  .  (2-18) 

Furthermore,  he  noted  that  the  outer  layer  has  a  wakelike  structure 

* 

and  its  eddy  viscosity  (p^)  scaled  with  6  and  U.  Thus, 

Ut  -  KPU6*  ,  (2-19) 

where  K  is  a  constant. 

Mellor  and  Gibson  [23]  combined  equations  (2-13)  and  (2-19)  into 
a  theory  of  equilibrium  wake  profiles,  which  is  in  good  agreement 
with  experimental  data,  as  shown  in  Fig.  6-7,  p.  479  of  Ref.  [19], 

The  equilibrium  concept  is  valid  not  only  for  positive  values  of  6  but 
also  for  negative  B,  as  shown  by  Stratford  [24],  and  Herring  and 
Norbury  [25],  for  example. 

Coles  [26]  combined  the  wall  component  and  the  wake  component 
with  the  wake  function  (W)  and  the  wake  parameter  (II)  in  order  to 
express  the  velocity  profile  in  the  form 


1 

K. 
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(2-20) 


The  wake  function  is  nor.'r lized  such  that  it  is  zero  at  the  wall  and 
attains  a  value  of  2.0  at  y  =  6.  As  shown  by  Coles,  the  wake  function 


- - — 
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is  approximated  by 

W(J)  *  2  sin2(||)  ,  (2-21) 

which  is  quite  av  curate  for  routine  use  as  shown  in  Fig.  6-8,  p.  481 
of  Ref.  [19].  The  wake  parameter  derived  by  Coles  and  Hirst  [27]  is 
approximated  by  the  empirical  expression 

n  «  0.8  (8  +  0.5)0,75  .  (2-22) 

This  formula  is  arranged  to  fit  the  theoretical  requirement  of  Mellor 
and  Gibson  [23]  i.e.,  the  wake  must  vanish  at  8  =  -0.5,  corresponding 
to  an  asymptotically  large  favorable  gradient.  For  a  very  large  8, 
where  there  are  no  equilibrium  data,  the  wake  parameter  was  approximated 
by  the  relation 


n  a  1  +  2.1  . 


(2-23) 
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CHAPTER  III 

PRELIMINARY  CONSIDERATIONS 


3.1  Definition  of  Cavitation  Number 

Cavitation  is  defined  as  the  vapor  and  gas-filled  voids  occurring 
in  a  liquid  when  the  local  pressure  falls  below  the  vapor  pressure, 
based  on  the  bulk  temperature  of  the  fluid.  The  Euler  number  which 
describes  this  state  may  be  written  as 


a  = 


p  -  p 


1/2  V, 


2  ’ 


(3-1) 


where  O  is  defined  as  the  cavitation  number,  P  and  U  are  a  reference 
pressure  and  velocity,  respectively.  The  vapor  pressure  (P  )  and 
liquid  mass  density  ( p)  correspond  to  the  bulk  temperature,  of  the 
liquid. 

The  cavitation  number  is  a  valuable  parameter  in  the  correlation 
of  cavitation  data. 

If  0  is  sufficiently  large,  there  is  no  evidence  of  cavitation; 
whereas  if  CT  is  sufficiently  small,  the  cavitation  can  be  well 
developed.  Limited  cavitation  occurs  at  an  intermediate  value  of 
the  cavitation  number  where  the  amount  of  cavitation  is  minimized. 

The  cavitation  number  corresponding  to  this  state  is  the  limited 
cavitation  number  (a^)  given  by 


a 
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i/2  P  U 
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(3-2) 
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3.2  Incipient  Cavitation  and  Desinent  Cavitation 

There  are  two  ways  to  determine  the  cavitation  number  at  a  constant 
velocity  in  a  water  tunnel.  The  first  method  is  to  establish  a 
velocity  in  the  test  section  at  a  noncavitating  condition  and  then 
decrease  the  tunnel  pressure  until  cavitation  occurs  on  the  test 
element.  The  cavitation  number  thus  determined  is  called  the  incipient 
cavitation  number  (a^).  The  second  method  is  to  establish  a  velocity 
in  the  test  section  at  a  condition  of  well-developed  cavitation  on 
the  test  element.  The  pressure  is  increased  until  the  cavitation 
disappears.  This  method  determines  the  desinent  cavitation  number  (o^) • 
Experience  shows  that  cr^  is  often  more  consistent  than  and  tends 
to  be  the  upper  bound  of  Therefore,  in  many  cavitation  experiments, 

cr^  is  adopted  as  the  limited  cavitation  number  (a^).  Differences 
between  incipient  and  desinent  cavitation  arc  due  to  cavitation 
hysteresis  [28]. 

3 . 3  Analysis  of  Isolated  Roughness 

Isolated  roughness  may  be  divided  into  two  categories,  namely, 
streamlined  or  non-separating  surface  protrusions  and  separating  such 
as  a  triangular  element.  For  the  former  cases,  cavitation  occurs  on 
the  surface  of  the  protrusion.  For  the  latter  case,  cavitation  occurs 
in  the  shear  layer  downstream  from  the  point  of  separation. 

The  frozen-streamline  theory  (FST)  developed  by  Holl  [29,30]  has 
been  used  to  analyze  the  major  characteristics  of  the  flow  about  a  two- 
dimensional  nonst  >araced  surface  protrusion.  Further  an  approximate 
method  developed  by  Holl  [4,29],  which  is  called  the  characteristic 
velocity  theory  (CVT),  was  able  to  correlate  the  limited  cavitation 
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number  more  successfully  than  FST,  for  a  non-streamlined  surface 
protrusion. 

3,4  Characteristic  Velocity  Theory 

Holl  [A]  developed  the  theory  based  on  the  assumption  of  frozen 
streamlines.  The  actual  flow  about  the  roughness  element  on  a  flat 
plate  In  rotational  flow  is  shown  in  Figure  3-1.  In  order  to  gain  a 
simple  expression  for  the  minimum  pressure  coefficient  of  the 
irregularity  in  rotational  flow  on  a  flat  plate  (Cp  )  an  approximate 
rotational  flow  is  constructed  and,  in  addition,  a  chl«c terist ic 
velocity  GO  is  introduced  as  shown  in  Fig.  3-2. 

Between  the  flow  at  infinity  along  the  wall  and  minimum  pressure 
point,  Bernoulli's  equation  is  valid  in  the  following  form: 


P 


P  .  _  +  U  1 

min  R  2  max 


:  can  be  expressed  in  the  form 

^min-o 


(3-3) 


C 

p 

min-o 


U  2 
,  max. 

U  „ 


U  2 

(— ) 
U_  • 


(3-4) 


If  we  assume  that  the  shapes  of  the  streamlines  are  the  same  in  both 
the  rotational  *ad  irrotational  flow,  then  ve  can  assume  that 
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U  2 
=  (,J5.ax 

VU’  ’ 
max 


(3-5) 


Substitution  of  equation  (3-5) 


into  equation  (3-4)  results  in 
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U  2 

c  -  (_£)  C’  ,  (3-6) 

D  ,  U  D  . 

tmin~o  00  min-o 


where  C'  is  the  minimum  pressure  coefficient  of  the  irregularity 

^min-o 

in  irrotational  flow  at  a  flat  plate.  If  cavitation  is  assumed  to 
occur  at  vapor  pressure,  i.e. ,  that  P  ,  _  =  P  ,  then  the  minimum 

pressure  coefficient  may  be  replaced  by  the  limited  cavitation  number 
on  a  flat  plate  and  equation  (3-6)  then  becomes 


°lfP  ~  a^P 


(3-7) 


where  is  the  limited  cavitation  number  for  the  roughness  in  the 

absence  of  a  boundary  layer,  that  is  h/6  =  cn.  Holl  then  assumed  that 
is  the  square  root  of  the  mean  value  of  the  square  of  the  velocity 


upstream  of  the  protrusion  over  some  characteristic  length,  L^.  Holl 
assumed  that  is  equal  to  h  so  that 


U 

c 


(3-8) 


For  turbulent  flat  plate  boundary  layers,  power  law  velocity 
profiles  are  reasonable  approximations  to  the  actual  velocity  profiles 
at  infinity  except  at  the  wall  where  the  slope  is  infinite.  Thu  we 
have 
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Employing  equation  (3-8)  in  equation  (3-7)  for  the  power  profiles 
yields 


£f  p 


1  ,h  “  a* 

H  (6} 


for  g<l 


(3-10) 


and 

°Up-  fl>r  l-i  -1  •  (3-n) 


The  shape  factor  (H)  defined  as  the  ratio  of  the  displacement  thickness 
(6*)  and  the  momentum  thickness  (0)  has  the  value 


H 


tn-4-2 

m 


(3-12) 


for  the  power  law  velocity  profiles  given  by  equation  (3-^). 

The  theory  may  be  employed  to  obtain  a  rough  estimate  of  a 
for  a  separated  protrusion  element. 
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3 . 5  General  Estimation  of  the  Limited  Cavitation  Number  on  a  Flat  Plate 

for  Isolated  Roughness 

We  would  expect  that  the  cavitation  number  for  an  isolated  roughness 
should  be  a  function  of  the  relative  roughness  (h/6)  and  the  boundary 
layer  shape  factor  (H)  ,  from  the  f rczen-st reamline  and  characteristic 
velocity  theories.  In  addition,  Reynold's  number  (Re^),  based  on  the 
height  of  roughness,  should  be  added  to  the  function,  according  to  the 
studies  by  Kermeen  and  Parkin  [31 J,  and  Arndt  [32].  Thus,  we  would 
expect  the  cavitation  number  for  isolated  roughness  on  a  flat  plate  to 
have  the  general  form 
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which  was  proposed  by  Holl  [29,30].  As  a  first  approach,  a  simple 
power  law  equation  for  an  isolated  roughness  on  a  flat  plate  may  be 
given  as 


„,h.a  ,US.b  1 

C£fp  ”  C(6)  (~\> )  d  » 

H 


(3-14) 


where  C,  a,  b,  and  d  are  constants  that  depend  on  the  shape  of  the 
irregularity. 


3 . 6  Superposition  Equation 

Consider  a  surface  irregularity  located  at  some  point  on  a  parent 

body  as  shown  in  Fig.  3-3.  The  minimum  pressure  (P  ,  „)  occurs  on 

7  r  min-R 

the  surface  irregularity  or  in  the  wake  of  the  irregularity  if  the  flow 

is  locally  separating  due  to  the  roughness.  The  pressure  difference 

P  .  -  P  can  be  written  as 

min-R  °° 


P  -  P 

min-R  ® 


P  -  P 


+  P  -  P 

min-R 


(3-15) 


where  P  is  the  pressure  in  the  vicinity  of  the  roughness  which  may  be 

thought  of  as  the  "local  pressure  at  infinity"  for  the  roughness  element. 

2 

Dividing  equation  (3-15)  by  1/2  pU  yields 


p  .  -  p  p  -  p  p  ,  -  p  „  : 

min-R _ “  _ _ mln-R  ,  U  . 

">  2  +  2  MI  ' 

1/2  p'J~  1/2  pU  1/2  pU 


(3-16) 


Figure  3-3.  A  Surface  Irregularity  on  a  Parent  Body. 
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where  U  is  the  velocity  at  the  edge  of  boundary  layer  in  the  vicinity 
of  roughness  which  may  be  thought  of  as  the  "local  velocity  at  infinity" 
for  the  surface  irregularity.  In  terms  of  pressure  coefficients, 
equation  (3-16)  can  be  written  in  the  following  form 


C 

Pmin-R 


C  + 
ps 


(1-C  )C 

pS  Pmin-o 


(3-17) 


Furthermore,  if  cavitation  is  assumed  to  occur  at  vapor  pressure, 
then  the  minimum  pressure  coefficients  may  be  replaced  by  cavitation 
numbers  and  equation  (3-17)  then  becomes 
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hr 


c  + 
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(l-c  )  a 

ps  Jofp 


(3-18) 


This  is  the  so-called  superposition  equation  and  was  derived  by 
Holl  [29]. 


The  equation  is  very  useful  in  estimating  the  cavitation  performance 
of  a  hydronautical  body  with  surface  irregularities,  provided  the  smooth- 
body  pressure  coefficients  (C  )  anc*  t^ie  limited  cavitation  number  of 
the  surface  irregularity  on  a  flat  plate  are  known. 

This  equation  has  been  revised  by  Amdt ,  Holl,  Bohn,  and  Bechtel 
[33]  to  include  possible  bubble  dynamic  effects.  The  revised  equation 


is 
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IR 


C  -  f;  +  (1-C  ) 

ps  ps 


(3-19) 


The  function  .Q  accounts  for  possible  bubble  dynamic  effects  and  is 
defined  as 


>  r~i  V > nt  -  ~Vl  j,  a. 


20 


Q, 


P  -  P  , 
v  min-r 

1/2  pl£, 


(3-20) 


where  P  ^  is  the  minimum  pressure  produced  by  the  roughness. 
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CHAPTER  IV 

EXPERIMENTAL  EQUIPMENT  AND  PROCEDURES 

4 . 1  General  Description  of  Water  Tunnel 

The  experiments  in  this  study  were  performed  in  the  Garfield  Thomas 
Water  Tunnel  Building  of  the  Applied  Research  Laboratory  located  at 
The  Pennsylvania  State  University.  The  12-inch  diameter  water  tunnel, 
facility  was  designed  to  accommodate  circular  (12- inch)  and  rectangular 
(4.5-inch  x  20-inch)  test  sections.  A  general  view  of  this  facility  is 
shown  in  Fig.  4-1.  The  rectangular  test  section,  shown  in  Fig.  4-2 
and  Fig.  4-3,  was  used  throughout  this  study. 

The  flow  velocity  can  be  changed  from  10  to  70  fps  by  a  variable 
speed  mixed-flow  pump.  The  pressure  in  the  test  section  can  be  varied 
from  3  to  60  psia  by  means  of  a  partially  filled  water  tank  having  an 
air  chamber  on  top  and  connected  to  the  tunnel  circuit  at  the  bottom. 
Alterations  of  the  air  pressure  in  the  tank  changes  the  pressure  in 
the  entire  system. 

The  gas  content  was  controlled  by  using  a  Cochrane  Cold-Water 
Degasifier  located  in  a  bypass  circuit.  The  .system  can  decrease  the 
total  gas  content  to  2  moles  of  air  per  10^  moles  of  water  (PPM) . 

For  more  detailed  information,  see  reference  [34]. 

4 . 2  Details  of  the  Test  Corif  igurat  ions 

The  experimental  arrangement  in  the  rectangular  test  section  was 
composed  of  a  triangular  protrusion,  the  base  for  the  protrusion  and  a 
curved-wall  element.  The  test  section  fitted  with  two  10- inch  by  20- inch 
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plexiglass  windows  with  an  additional  2- inch  by  20- inch  window  on 
the  bottom  panel.  This  arrangement  is  sketched  in  Fig.  4-4. 

4.2.1  Curved  Wall.  The  curved  wall  was  made  from  plexiglass. 

The  dimensions  of  the  wall  are  shown  in  Fig.  4-5.  In  order  to  measure 
the  velocity  profiles  at  each  location  on  the  flat  plate,  seven  holes 
were  drilled  in  the  wall  for  a  total  head  tube.  These  locations  are 
designated  STA. I,  II,  III,  IV,  V,  VI  and  VII,  as  shown  in  Fig.  4-6. 
Additional  details  of  the  design  are  presented  in  Appendix  A. 

4.2.'1  Triangular  Protrusion.  The  triangular  protrusions  used 
in  this  study  were  made  from  plexiglass.  A  photograph  of  the  pro¬ 
trusion  is  shown  in  Fig.  4-7.  The  protrusions  were  of  the  same  shape 
as  those  used  by  Holl  [29].  Eight  roughness  heights  were  employed  in 
the  tests;  namely,  1/100,  2/100,  3/100,  5/100,  7/100,  10/100,  15/100 
and  30/100  inches  and  were  mounted  on  the  base  one  at  a  time  for  each 
experiment.  The  roughness  locations  correspond  to  the  position  of 
the  leading  edge  of  the  total  head  tube.  Additional  details  of  the 
design  are  presented  in  Appendix  A. 

4 . 3  Meas  urement  of  the  Velocity  :i  t  Infinity  in  the  Test  Section 

Infinity  in  the  test  section  was  uefined  as  6  40  inches  upstream 
from  the  leading  edge  of  the  curved  wall  for  all  experiments.  The 
velocity  was  determined  by  the  difference  between  the  static  pressure 
at  the  infinity  point  with  the  curved  wall  (P^)  and  the  static  pressure 
far  upstream  of  the  infinity  point  (P  ) .  The  relationship  between  the 

static  pressure  difference  (P  -  P'l  and  the  dynamic  pressure  at 

o  °°  ■ 

1  7 ' 

infinity  with  the  curved  wall  (—  PU  )  was  measured  three  times  by 
using  a  pitot  tube.  However,  since  cavitation  occurred  at  the  pitot 
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tube  for  flow  velocities  greater  than  40  fps,  the  measurements  were 
obtained  for  velocities  from  20  to  38  fps.  The  sketch  of  the  measure¬ 
ment  stations  is  shown  in  Fig.  4-8. 

4 . 4  Measurement  of  the  Pressure  Coefficient 

The  pressure  coefficients  for  the  test  section  were  determined 
by  using  a  total  head  tube  and  the  static  taps  on  the  flat  plate.  The 
static  taps  are  located  between  the  measuring  stations  on  a  smooth 
plexiglass  plate,  as  shown  in  Fig.  4-9.  The  measurements  by  the  total 
head  tube  and  the  static  taps  were  not  made  at  the  same  time  because 
the  flow  is  disturbed  by  the  presence  of  a  total  head  tube.  First 
the  static  pressure  at  each  station  was  measured  for  different  flow 
velocities  ranging  from  20  fps  to  60  fps  and  then  the  dynamic  pressure 
at  each  station  was  measured. 

4 . 5  Measurement  of  the  Velocity  Profile 

Whenever  the  measurements  of  the  mean  velocity  profile  for  each 
station  were  made,  the  positions  at  which  the  tip  of  the  total  head  tube 
just  touches  the  plate  for  both  inward  and  outward  movement  of  the 
traversers  were  carefully  determined  with  the  use  of  a  microscope. 

This  was  essential  in  order  to  know  the  back-lash  of  the  traversing 
gear  and  also  to  determine  the  initial  measuring  point.  This  procedure 
was  of  major  importance  in  establishing  accurate  locations  of  the  total 
head  tube. 

The  velocity  profiles  were  measured  by  traversing  the  boundary 
layer  with  a  probe  which  was  fabricated  from  a  tube  of  rectangular 
cross  section  with  an  external  size  of  C.018  by  0.083  Inches.  The 
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traversing  device  pictured  in  Fig.  4-10  allowed  the  positioning  of 
the  total  head  tube  within  0.001  inch. 

For  each  velocity  profile,  the  total  head  tube  position  was 
zeroed  on  the  plate  and  then  moved  outward  in  increments  of  0.010  inch 
for  a  few  initial  measuring  points,  after  which  it  was  moved  in  steps 
of  either  0.020  or  0.050  inch  until  the  edge  of  the  boundary  layer 
was  reached,  in  other  words,  until  the  total  pressure  became  constant. 

The  measurements  were  repeated  as  the  tube  was  moved  inward  toward 
the  wall  in  the  same  increments  as  employed  in  the  outward  motion.  A 
third  data  set  was  obtained  as  the  probe  was  again  moved  outwards. 

The  final  profiles  were  determined  from  the  average  of  these  three 
sets  of  data.  While  making  measurements,  the  tunnel  pressure  was 
maintained  at  40  psia  and  the  water  temperature  was  recorded  before 
and  after  each  profile  run. 

4 . 6  Visual  Observation  of  Desinent  Cavitation 

Before  observing  cavitation,  the  gas  content  in  the  water  tunnel 
was  always  measured.  If  the  value  was  more  than  4  p.p.m. ,  purging 
the  water  of  gas  was  done  before  the  experiment  by  using  the  tunnel 
degasser.  The  gas  content  was  measured  by  a  "Thomas"  Van  Slyke 
Manonetric  Apparatus  using  a  10  cc  sample. 

As  previously  mentioned,  there  are  two  methods  by  which  the  critical 
cavitation  number  may  be  determined,  i.e.,  inception  and  desinence. 
However,  the  desinence  method  was  used  to  determine  the  limited  cavitation 
number  during  the  entire  test  because  of  repeatability. 
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At  the  start  of  the  observation,  a  strobe  light  was  used  to  aid  in 
visually  observing  cavitation.  When  cavitation  occurred  along  the  top 
edge  of  the  triangular  protrusion,  which  appeared  as  luminous  streaks 
originating  several  inches  downstream  of  the  protrusion,  another  type 
of  cavitation  appeared  at  both  end3  of  the  protrusion.  These  two  types 
of  cavitation  joined  together  downstream.  Special  care  was  taken  to 
see  that  the  cavitation  caused  by  the  ends  was  not  mistaken  for  that 
caused  by  the  top  edge  of  the.  protrusion. 

Cavitation  disappearance  was  defined  as  the  absence  of  cavitation 
at  the  top  edge  of  the  protrusion,  for  less  than  39  seconds.  Throughout 
the  experiment,  the  observations  were  done  by  the  same  person.  The 
tunnel  pressure  was  recorded  at  the  defined  infinity  point.  For  each 
velocity  setting,  this  pressure  was  lowered  and  raised  seven  times, 
with  the  final  deslnent  cavitation  number  being  based  on  the  average 
pressure  obtained  by  this  procedure.  In  addition  to  the  pressure 
measurements,  the  gas  content  and  temperature  were  measured  before  and 
after  each  run.  Since  one  side  of  the  transducer  was  open  to  the 
atmosphere  while  taking  the  pressure  measurements,  it  was  necessary 
to  record  the  barometric  pressure  several  times  each  day. 

A  photograph  of  the  cavitation  phenomenon  which  occurred  on  the 
triangular  protrusion  is  shown  in  Fig.  4-11. 
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CHAPTER  V 

PRESENTATION  AND  DISCUSSION  OF  RESULTS 

5. 1  Velocity  at  Infinity  In  the  Test  Section 

As  indicated  previously,  infinity  in  the  test  section  was  6.40 

inches  upstream  from  the  leading  edge  of  curved  wall  (see  section  4.3). 

The  velocity  at  infinity  cannot  be  measured  directly  by  using  the 

reference  static  pressure  in  the  test  section  because  of  the  presence 

of  the  curved  wall  in  the  test  section. 

The  relationship  between  the  static  pressure  difference  (P  -  P^) 

2 

and  the  dynamic  pressure  at  infinity  (1/2  pU^  )  was  empirically 
established  from  a  set  of  experiments  (see  Fig.  4-8)  and  is  of  the 
form 

1/2  pu;2  =  0.9708  (Pq  -  P^),  (5-1) 

where  U'  and  P^  are  the  velocity  and  pressure  at  infinity  with  the 
curved  wall,  respectively,  as  sketched  in  Fig.  5.1.  The  values  of 
the  dynamic  pressure  and  static  pressure  difference  are  tabulated  in 
Table  5.1. 

The  accuracy  of  the  coefficient  (0.9708)  was  investigated  by 
applying  Bernoulli's  equation  and  also  by  utilizing  the  Douglas- 
Nsumann  computer  program.  Tne  relationship  between  the  dynamic 
pressure  and  static  pres  jure  without  the  curved  wall  can  be  established 
from  Bernoulli's  equation  by  knowing  the  ratio  of  cross-sect ional  area 
at  the  infinity  point  to  that  far  upstream  of  the  infinity  point.  Thus 


TABLE  5-1 


*  2 

RELATIONSHIP  BETWEEN  1/2 pU^  AND  (P^P*) 


DATA 

SET  1 


1/2  pi/ 2 
(volts) 

P 

o  00 
(volts) 

l/2pt/2 
(volt  s) 

P  -K 
o  “ 

(vo Its) 

1/2 PUp 

(vo)  ts) 

P  -P’ 
o  00 

(volts) 

Velocity 

(fps) 

20 

0.9446 

0.9850 

0.9255 

0.9780 

- 

— 

23 

1.2435 

1.3165 

1.2515 

1.3230 

1.2616 

1.3610 

25 

1,4859 

1.5537 

1.5073 

1.5715 

1.5596 

1.6580 

28 

2.0154 

2.0850 

2.0170 

2.0912 

1.9721 

2.0800 

30 

2.3290 

2.4206 

2.3027 

2.4085 

2.3467 

2  4400 

33 

2.7580 

2.8450 

2.7495 

2.8650 

2.6712 

2.8113 

35 

3.0054 

3.1231 

3.0113 

3.1460 

2.9897 

3.1320 

38 

3.4888 

3.6240 

3.5039 

3.6250 

3.4886 

3.6300 

l/2pt/2/P 

-Pr  0.9663 

o  00 

0. 

9693 

0. 

,9764 

Average  value  of  1/2  PUp/P^P^  -  0.9708* 


*Standard  deviation  *  0,0052 
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equation  (5-5)  into  equation  (5-6)  yields 

p  -  Pro 

CP  “  n004lTp~'-“O 

o  00 


(5-7) 


In  order  to  check  the  reliability  of  equation  (5-7),  the  pressure 
coefficient  was  also  calculated  from  the  velocity  at  the  edge  of  the 
boundary  layer  (U) ,  shown  in  Fig.  5-2,  which  was  measured  by  the  total 
head  tube.  The  pressure  coefficient  is  expressed  in  the  form  of 

U  ^ 

c  =  1  -  (— -)  .  (5-8) 

P  u® 

As  shown  in  Fig.  5-3,  these  pressure  coefficients  derived  from  equations 
(5-7)  and  (5-8)  are  in  fairly  good  agreement.  So  the  pressure  coefficient 
obtaineu  by  using  the  static,  tap  was  used  during  the  entire  experiment 
and  tabulated  in  Table  5-2. 

However,  in  comparison  with  the  calculated  pressure  coefficient 
which  was  used  for  designing  the  carved  wall,  the  measured  coefficients 
differ  significantly.  The  difference  between  the  experimental  and 
calculated  values  is  greater  for  the  pressure  gradient  than  for  the 
pressure  distribution.  These  two  pressure  gradients  are  shown  in 
Fig.  5-4.  Since  STA.  Ill  and  V  have  almost  the  same  pressure  gradient 
as  that  of  STA.  IV,  the  experiment  for  observing  cavitation  was  not 
carried  out  for  these  two  stations. 

5 . 3  Boundary  Layer  Parameters  on  the  ^late 

5.3.1  Boundary  Layer  Thickness.  The  boundary  layer  thickness 
was  defined  co  be  the  distance  from  the  wall  where  the  mean  local 
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velocity  is  99%  of  the  velocity  at  the  edge  of  the  boundary  layer. 

The  boundary  layer  thickness  (5)  for  each  station  was  obtained  from 
the  velocity  profiles.  Furthermore,  these  values  were  corrected  with 
respect  to  the  velocity  at  infinity  (U^)  equal  to  40  fps  at  a  temperature 
of  74°F. 

The  equation  used  for  the  correction  is 


6  ■ 


0. 37x 
(Ux/v)1/5 


(5-8) 


which  describes  the  growth  of  a  turbulent  boundary  layer  thickness  on  a 

flat  plate.  The  corrected  values  are  shown  in  Table  5-3  and  Fig.  5-5. 

5.3.2  Momentum  Thickness.  Other  parameters  such  as  the  displace- 
* 

ment  thickness  (o  )  and  the  momentum  thickness  (9)  are  given  by 


(5-9) 


and 


(i  -  x> 


(5-10) 


respectively,  and  these  could  not  be  obtained  from  the  velocity  profiles. 

The  minimum  ratios  between  the  mean  velocity  (u)  and  the  velocity  at 

the  edge  of  boundary  layer  should  at  least  be  less  than  0.6  for  this 

experiment  in  order  to  obtain  the  momentum  thickness  (B)  and  the 

* 

displacement  thickness  (6  ).  Unfortunately,  the  minimum  values 
obtained  by  the  experiment  are  more  than  0.8.  To  circumvent  this 
problem,  the  method  of  E.  Truckenbrodt  was  employed  to  obtain  the 

■k 

momentum  thickness  (6)  and  the  displacement  thickness  (6  ).  in  so  doing, 
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TABLE  5-3 

BOUNDARY  LAYER  THICKNESS  AT  EACH  STATION 

6+f  at  Uoo=,40fps 


STA 

X 

( inches) 

U 

(fps) 

Temp . 
(°F) 

S+ 

( inches) 

T>74°F 
( inches) 

I 

1.7 

43.07 

76.5 

0.125 

0.124 

II 

4.3 

45.22 

69.5 

0.154 

0.156 

III 

5.6 

42.91 

75.0 

0.160 

0.162 

IV 

7.0 

42 . 60 

74.0 

0.165 

0.167 

V 

8.7 

42.52 

74.0 

0.180 

0.183 

VI 

11.1 

41.89 

74.0 

0.230 

0.233 

VII 

13.6 

40.84 

75.0 

0.370 

0.374 

+  6 

as  measured  for 

each  TJ. 

++  6  at  U«j*40fps,  T*74°F  was  obtained  by  using  equation  (5-8)  and 

Fig.  5-2  with  respect  to  the  velocity  at  infinity  (Um)  equal 
to  40  fps  at  the  temperature  of  74°f. 
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two  computer  programs  based  on  the  method  of  E.  Truckenbrode  [35]  were 
used.  One  is  the  MLB  program  which  was  used  mainly  for  calculating 
the  above  parameters,  and  the  other  is  an  auxiliary  for  the  MLB  program, 
the  RJD  program. 

First  of  all,  the  MLB  computer  program  was  used  to  calculate  the 
momentum  thickness  (0)  and  the  boundary  layer  thickness  (6).  The 
origin  of  the  boundary  layer  was  shifted  back  and  forth  in  order  to 
obtain  the  best  agreement  between  calculated  and  experimental  values 
of  the  boundary  layer  thickness.  After  these  adjustments,  the  best 
agreement  between  calculated  and  measured  values  of  the  boundary 
layer  thickness  was  obtained  when  the  origin  for  the  calculated 
values  was  one  half  chord  length  upstream  of  the  leading  edge  of  the 
curved  wall,  as  shown  in  Fig.  5-5.  Then  the  momentum  tnicknesis  (.6) 
was  calculated  for  the  boundary  layer  thickness  employing  the  origin 
established  by  the  above  method. 

5.3.3  Shape  Factors.  Since  the  displacement  thickness  (6  ) 
could  not  be  obtained  by  the  MLB  computer  program,  the  RJD  computer 
program  was  used  to  calculate  the  displacement  thickness  and  shape 
factor  (H) .  By  knowing  the  pressure  distribution,  the  momentum 
thickness  and  a  new  shape  factor  (L)  can  be  calculated  by  using  the  RJD 
computer  program.  The  shape  factor  (1)  is  expressed  In  the  form 

f  dH  ,  (5-11) 

L  * - ~ 

J(H  -  1) H 


where 
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H 


and  the  energy  thickness  is  given  by 


** 


-2  - 

j_,  H. 

/  u 


(1  -  T-)  w  dy 

n 


(5-12) 


Initially,  the  momentum  thickness  was  calculated  for  the  origin 
established  for  the  boundary  layer  thickness.  As  can  be  seen  from 
Fig.  5-6,  the  difference  between  the  momentum  thickness  calculated  by 
the  two  computer  programs  is  not  insignificant.  This  is  due  to  the 
fact  that  the  MLB  program  utilizes  a  somewhat  different  skin  friction 
coefficient  (C.)  than  that  employed  in  the  RJD  program  which  is  based 
on  used  by  Truckenbrodt .  However,  the  calculation  of  the  shape 
factor  (H)  via  the  determination  of  the  other  shape  factor  (L)  shown 
in  Fig.  5-7  was  not  sensitive  to  these  differences  in  the  momentum 
thickness  (8)  as  evidenced  by  the  nearly  equal  values  of  H  tabulated 
in  Table  5-4.  This  lack  of  sensitivity  to  variations  in  H  may  be  due 
to  the  fact  that  the  slope  of  G  versus  X  shown  in  Fig.  5-6  is  nearly 
the  same  for  both  methods  of  calculation. 

It  is  felt  that  the  determination  of  6  and  9  employing  the  method 
of  E.  Truckenbrodt  as  described  above  is  fairly  accurate  oince  it  is 
based  upon  an  origin  determined  by  the  experimental  values  of  <5  as 
described  in  Section  5.3.2  (see  Table  5-5). 

5.3.4  Shear  Stress  at  the  Wall.  The  experimental  values  for 

the  shear  stress  at  the  wall  (t  )  can  be  approximated  by  the 

w 

empirical  formula 


mm 


Figure  5-6.  Momentum  Thickness  in  the  Test  Section 


HI  -B'Hllh--- 


z .zr  \ 


—  TRUCKENBRODT 

—  GARNER 

—  GRUSCHWITZ 


Figure  5-7,  Relation  Between  the  Shape  Factors 
(H)  and  (L)  for  Turbulent  Flow. 
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TABLE  5-5 

BOUNDARY  LAYER  PARAMETERS  AT  EACH  STATION  * 


STA 

X 

c 

6 

(  inches) 

6* 

( inches) 

0 

( inches) 

H 

I 

0.10 

0.124 

0.0122 

0.0091 

1.34 

II 

0.25 

0.156 

0.0148 

0.0111 

1.33 

III 

0.33 

0.162 

0.0162 

0.0123 

1.32 

IV 

0.41 

0.167 

0.0183 

0.0137 

1.335 

V 

0.51 

0.183 

0.0213 

0.0159 

1.34 

VI 

0.65 

0.233 

0.0270 

0.0197 

1.37 

VII 

0.80 

0.374 

0.0354 

0.0255 

1.39 

*These 
40  fps 

parameters  were  calculated 
(see  Section  5.3.1) 

for  a  velocity 

at  infinity 

(UJ  of 

-«*»*■ 


fla. 


which  was  given  by  Ludwieg- Tillman  [36] „  The  shear  stress  at  the 
wall  was  calculated  by  substituting  the  boundary  layer  parameters  which 
were  obtained  in  Section  5.3.2  and  5.3.3  into  equation  (5-13).  Further¬ 
more,  Clauser's  equilibrium  parameter  ( 0) ,  the  shape  factor  (G)  and 
the  wake  parameter  (II)  were  calculated  from  equations  (21-14)  ,  (2-18) 
and  (2-22),  respectively,  and  are  tabulated  in  Table  5-6  and  Table 
5-7. 

5.4  Cavitstion  Number 

5.4.1  Cavitation  Number  for  the  Case  of  Zero  Pressure  Gradient. 

The  cavitation  number  (a^  )  for  the  case  of  zero  pressure  gradient 
was  measured  at  STA.  IV,  and  is  tabulated  in  Table  5-8.  The  effect  of 
the  relative  height  (h/ 6)  has  been  known  co  be  the  most  dominant  factor 
for  this  case  and  the  relationship  between  the  cavitation  number  (a^  ) 

and  the  relative  height  for  this  investigation  (h/6)  is  shown  by  the 
experimental  data  in  Fig.  5-8.  A  comparison  of  the  present  result 
with  that  of  Ref.  [29]  was  carried  out  and  is  shown  in  Fig.  5-°. 

For  these  two  investigations,  Fig.  5-8  and  Fig.  5-9  show  that  the 
relationship  between  the  log  of  the  cavitation  number  and  the  log 
of  the  relative  height  is  .Linear.  By  neglecting  the  effect  of  Reynold's 
cumber,  these  relationships  can  be  expressed  in  the  form 

/a.  a 

o,  =  constant  (— ) 

*- 1  p  ^ 


(5-14) 


56 


TABLE  5-6 


WALL  SHEAR  STRESS 

AT  EACH  STATION 

STA 

X 

U 

6 

H 

T  * 
w 

c 

(ft/sec) 

(laches) 

(Ib/ft 

I 

0.1C 

43.30 

0.0091 

1.34 

3.366 

II 

0.25 

45.31 

0.0111 

1.33 

3.508 

III 

0.33 

45.83 

0.0123 

1.32 

3.537 

IV 

0.41 

46.04 

0.0137 

1.335 

3.381 

V 

0.51 

45.92 

0.0159 

1.34 

3.209 

VI 

0.65 

45.17 

0.0197 

1.37 

2.811 

Vi  I 

0.80 

43.69 

0.0255 

1.39 

2.401 

*  These 

values 

were  calculated  by 

using  equation  (5- 

-13)  (see 

Section 

5.3.4). 
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TABLE  5-8 


LIMITED 

CAVITATION 

DATA  STA. 

IV 

dx 

*  0) 

h 

U 

Gas 

Content 

Water 

Temp. 

Standard 
Deviat ion 

6 

Cfps) 

(fps) 

(ppm) 

(°F) 

%>* 

0f  °tfp 

0.058 

29.77 

34.27 

2.90 

0.270 

0.004 

0.060 

38.82 

44.69 

to 

68.  C 

0.292 

0.017 

0.063 

60.36 

69.48 

3.28 

0.322 

0.013 

0.290 

29.88 

34.39 

1.50 

0.646 

0.011 

0.299 

38.59 

44.42 

to 

70.0 

0.666 

0.015 

0.312 

61.26 

70.51 

2.93 

0.676 

0.004 

0.531 

29.73 

34.22 

1.98 

0.912 

0.026 

0.599 

38.38 

44.18 

to 

67.5 

0.901 

0.030 

0.625 

60.27 

69.37 

3.12 

0.872 

0.014 

1.736 

29.46 

33.91 

2.13 

1.836 

0.090 

1.790 

38.06 

43.81 

to 

68.0 

1.84  7 

0.144 

- 

- 

- 

4.50 

- 

- 

*  These  are  mean  values  for  seven  measurements 
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Figure  5-8.  Variation  of  Cavitation  Number  with  Relative  Height. 
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The  power  number  '  a‘  and  the  correlation  coefficient  (R)  were  calculated 


by  using 


n  E  iogCJ  E  log(ffa  ) 

t-  ,  ,hx  x  i»l  i  i=l  v  i 

8<«  1  8  '  Wp\ - o - 


"  FiogA)  - 
1-lL  S-  1 


2  i*l  1 


(5-15) 


3  Clog(|) 

°log(a^p) 


(5-16) 


where  a  was  obtained  by  the  least  squares  method,  a.  -h.  and  o,  ,  . 

h  lo*(i>  1°8(0wP) 

are  standard  deviations  of  log(-g)  and  log(a^p)  ,  respectively. 

'a'  is  0.53  for  the  correlation  coefficient  0.9969  for  the  present 
experiment  whereas  'a'  has  the  value  0.40  for  the  experiment  of  Ref. 

[29].  This  difference  in  *a'  for  the  two  experiments  may  be  due  to  gas 
content  effects.  The  present  experiments  were  carried  out  with  the 
gas  contents  between  1.91  ppm  and  5.52  ppm.  On  the  other  hand,  the 
experiments  of  Ref.  [29]  were  conducted  for  gas  contents  of  3.46  ppm 
to  15.66  ppm.  For  comparison  purposes  the  data  of  Ref.  [29]  were 
corrected  for  different  gas  contents.  The  equation  used  for  the 


correction  is 


K(aR  -  3.71) 


at>  =  a_  - 

R-corr  R 


1/2  PU 


(5-17) 
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where  k  is  0.104,  a  is  a  limited  cavitation  number  of  Ref.  [29), 

R 

and  a  is  a  gas  content  of  Ref.  [29],  The  derivation  of  equation  (5-17) 
is  demonstrated  in  Appendix  B»  The  corrected  cavitation  numbers  are 
shown  in  Fig.  5-10,  which  demonstrates  that  the.  slope  of  the  variable 
between  the  log  of  the  corrected  cavitation  number  to  the  log  of  the 
relative  height  becomes  almost  the  same  as  the  present  one.  It 
might  be  said  that  0.53  is  a  more  universal  value  for  a  specific  gas 
content,  l.e.  »  ot  =*  3.71  ppm  than  the  value  of  0.40  employed  in  Ref. 

[33  ] .  Thus 


constant 


0.53 


(5-18) 


However,  the  effect:  of  Reynold's  number  based  on  boundary  layer 
thickness  should  be  accounted  for  to  obtain  a  more  accurate  relation¬ 
ship.  This  will  be  discussed  in  Section  5.4.3. 

5 . 4 . 2  Effect  of  Pressure  Gradient  on  Cavitation.  The  cavit  at  ion 
numbers  for  the  case  ol  favorable  and  unfavorable  pressure  gradient 
were  measured  at  STA.  I.  and  STA.  VII  and  are  tabulated  in  Tables  5-9 
and  5-10,  respectively.  Fig.  5-11  shows  that  the  relationships 
between  the  log  of  the  cavitation  number  and  the  log  of  the  relative 
height  are  almost  linear.  The  slopes  are  0.576  for  the  correlation 
coefficient  (R)  of  0.9885  for  the  case  of  a  favorable  pressure  gradient, 
and  0.481  for  R  equal  to  0.9882  for  the  case  of  an  unfavorable  pressure 
gradient.  The  cavitation  data  for  zero  pressure  gradient  lies  in 
between  the  data  shown  in  Fig.  5-11  and  never  exceeds  those  for  the 
case  of  unfavorable  pressure  gradient.  In  other  words,  the  cavitation 
numbers  tend  to  increase  with  the  pressure  gradient. 
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TABLE  5-9 


LIMITED 

CAVITATION  DATA, 

STA.  I  <-£■  <  0) 
dx 

U<x, 

Gas 

Water 

h 

U 

Content 

Temp.  a  * 

(Op)  ^fp 

6 

(fps) 

(fps) 

(ppm) 

0.079 

30.07 

35.53 

2.60 

0.247 

0.081 

38.82 

42.07 

to 

68.0 

0.265 

0.085 

60.89 

65.98 

3.90 

0.310 

0.156 

29.60 

32.07 

2.91 

0.503 

0.161 

37.49 

40.62 

to 

67.0 

0.545 

0.168 

59. 39 

64.35 

4.05 

0.603 

0.548 

29.31 

31.76 

1.98 

0.886 

0.565 

37.89 

41.06 

to 

68.5 

0.886 

0.590 

59.01 

63.88 

2.72 

0.946 

0.782 

29.63 

32.10 

3.85 

1.040 

0.806 

38.31 

41.51 

to 

69.5 

1.065 

0.841 

59.32 

64.27 

4.50 

1.107 

2.347 

29.41 

31.87 

2.52 

2.130 

2.420 

37.79 

40.95 

to 

71.5 

2.125 

2.526 

47.98 

51.99 

3.52 

2.470 

Stcindard 
Devi  at:  ion 

°f  °*fp 


0.016 

0.005 

0.005 


0.010 

0.002 

0.006 


0.017 

0.006 

0.007 


0.014 

0.013 

0.025 


0.061 

0.053 

0.056 


*  These  are  mean  values  for  seven  measurements. 
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TABLE  5-10 


LIMITED  CAVITATION  BATA,  STA. 


VII  ( 


4r 

dx 


>  0) 


u 

00 

(fps) 

Gas 

Water 

Standard 

h 

6 

U 

(fps) 

Content 

(ppm) 

Temp. 

(°F) 

o  * 

2fp 

Devi at  ion 

o£  5«P 

0.026 

36.77 

42.35 

2 .25 

0.268 

0.006 

0.027 

49.30 

53 . 84 

to 

71.0 

0.281 

0.001 

0.02!? 

60.83 

66. 4  A 

3.15 

0.281 

0.009 

0.051 

30.07 

32.84 

3.30 

0.465 

0.008 

0.053 

38.82 

42.40 

to 

66.0 

0.506 

0.003 

0.055 

61.10 

66.  74 

5.21 

0.47/ 

0.010 

0.180 

79.82 

32.57 

3.98 

0.750 

0.011 

0. 186 

37„  85 

41.34 

to 

72.0 

0.787 

0.010 

0.194 

59.92 

65.45 

5.52 

0.851 

0.005 

0.389 

30.01 

32.78 

3.15 

1.012 

0.030 

0.401 

38.87 

42.45 

to 

70.0 

0.985 

0.012 

- 

— 

- 

5.32 

- 

- 

0.778 

29,46 

32.18 

2.36 

1.612 

0.082 

0.802 

37.89 

A1.53 

to 

70.0 

1.58G 

0.035 

0.837 

48.36 

32.82 

3.81 

1.697 

0.020 

*  These  are  mean  values  for  seven  measurements 


Variation  of  Cavitation  Number  with  Relative  Height  for  Finite  Pressure  Gradients 
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According  to  equation  (3-7)  used  in  the  characteristic  velocity 
theory,  that  is, 


V ’ <->  V' 


(3-7) 


the  velocity  profile  for  a  favorable  pressure  gradient,  as  shown  in 
Fig.  2-1,  i.  st  promote  the  onset  of  cavitation.  However,  this  investi¬ 
gation  showed  that  the  cavitation  number  increased  with  the  pressure 
gradient,  which  is  in  contradiction  with  the  characteristic  velocity 
theory.  This  discrepancy  may  be  due  to  the  complicated  flow,  i.e., 
separated  flow,  caused  by  the  triangular  protrusion. 

As  discussed  in  Section  2.3,  the  shape  factors  (H)  and  (G) ,  and 
wake  parameter  (II)  can  be  considered  as  factors  to  be  related  to 
pressure  gradient.  In  order  to  obtain  the  most  consistent  relationship 
between  cavitation  number  and  one  of  these  three  factors,  the  calculation 
was  done  by  using  th?  present  data  and  did  not  include  the  data  of 
Holl  [29]  .  Fig.  5-12  shows  that  power-law  correlations  of  <7^  with 
G  and  II  could  be  obtained  fur  the  present  data.  The  correlation 
coefficient  for  each  h/5  and  the  average  percent  error  in  slope  are 
0.993  and  18.0%,  respectively ,  tor  correlation  of  ^  with  G  and 

0.9901  and  21.4%  for  that  of  0„,  with  II.  It  is  seen  that  the 

x-fp 

correlation  of  0„  .  with  G  is  more  consistent  than  that  of  ft. 
tip 

However,  II  may  be  an  alternative  parameter  to  G  whereas  the  parameter 

H  may  be  primarily  limited  to  the  case  dp/dx  =  0. 

Thus,  the  relationship  between  the  cavitation  number  (o„  ) 

■  fp 

and  the  shape  factor  (G)  Is  expressed  for  the  same  relative  height 
and  Reynold's  number  in  the  form 


NOMINAL  VELOCITY  =  40  fps 


Figure  5-12.  Variation  of  Cavitation  Number  with  H,  II  and  G. 


69 


=  constant  (G) 


1.25 


(5-19) 


5.4.3  Effect  of  Reynold's  Number  on  Cavitation.  As  indicated 
in  Section  3.5,  Reynold's  number  based  on  the  boundary  layer  thickness 
i3  an  important  parameter.  However,  the  influence  of  Reynold's 
number  is  less  dominant  than  the  relative  height  or  the  shape  factor 
(G) ,  as  shown  in  Figs.  5-8  and  5-11.  The  relationship  between  the 
cavitation  number  and  the  Reynold's  number  for  constant  relative 
height  and  shape  factor  (G)  has  the  form 


,U6. 

:=  constant  ( — ) 
£fp  V 


0.194 


(5-20) 


Since  the  exponent  0.194  is  close  to  the  va" ue  0.20  determined  by 
Arndt  et  al.  [33]  ,  the  value  0  20  was  used  throughout  this  analysis. 
Hence 


°Hfp 


constant 


0.20 


(5-21) 


5.4.4  The  Cavitation  Law  for  Two  Dimensional  Triangular  Protrusions . 
The  cavitation  law  for  a  two  dimensional  triangular  protrusion  can 
be  expressed  in  the  form 


0.53 


°cr  "  C(T> 
~f  o  0 


A 

V 


0.2 


1.25 


(G) 


(5-22) 


When  applied  to  the  data  from  this  investigation  the  above  correlation 
gives  a  good  esc.i  tale  of  the  cavitation  number  with  an  average  error 
of  10. 37,  with  respect  to  the  measured  values.  The  equation  is  compared 
with  exper imenta1  data  in  Fig.  5-13. 
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i&'t 


As  shown  in  Fig.  5-14,  a  fair  agreement  wac  obtained  with  the 
results  of  Holl  for  -  0,  when  corrected  to  a  gas  content  of  3.71 
ppm  as  described  in  Appendix  B, 

5.5  Development  of  the  Characteristic  Velocity  Theory 


As  discussed  in  Section  3.4,  O ^  is  expressed  as 


3mp-s(!>"V  ,or 


°MP  * 1  - 1  ^  -  S>  °ifp  f°r  i  -  1 


(3-10) 


(3-11) 


where  ai,  is  the  limited  cavitation  number  for  the  rou;  hness  in  the 
xfp 

absence  of  a  boundary  layer,  that  is,  h/<S  =  For  the  two  dimensional 

triangular  roughness,  was  assumed  to  be  equal  to  that  of  a  sharp- 

edged  flat  plate  obtained  by  J.  W.  Holl  and  C.  B.  Baker  of  the  Applied 
Research  Laboratory  at  The  Pennsylvania  State  University,  and  is  shown 
in  Fig,  5-15.  By  the  method  of  least  squares,  ^  was  calculated  to 


O'  -  0.0214  Re, 
2.fp  h 


0.424 


(5-23) 


where  Re^  is  the  Reynold's  number  based  on  the  height  of  the  roughness. 

Examples  of  the  application  of  this  theory  based  on  equation 
(5-23)  are  shown  in  Fig.  5-8  for  the  nominal  velocities  40,  50  and 
bO  fps.  As  shown  in  Fig.  5-8,  it  is  seen  that  calculated  by 

the  characteristic  velocity  theory  is  in  fair  agreement  with  the 
experimental  data  for  0.0b  <  h/S  <  1.79,  which  was  the  investigated 
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HOLL  (1958)  DATA 

o  PHASE  I  FORWARD  STA. 
A  CENTER  STA. 

•  PHASE  n  FORWARD  STA. 

*  REAR  STA. 


I 


ai  fp 


h/6 

(U6/2-) 

__G_ 

0.429—3.400 

3.66x10^-8.67x10? 

6.52E 

0.139-1.100 

9.90x107-2.64x10"! 

6.47C 

C. 035-0. 524 

2.14x10^-6.81x107 

6.510 

0.014-0.213 

5.23x10 -1.51X105 

6.400 

CORRECTED  TO  GAS  CONTENT  OF  3.71  pom 
°lfp  =  0-0157  (h/M°- 53  (U6/j/)°-20  (G)h25 


ooa(^r  «c.iz5 


1000 


Figure  5-14.  Cavitation  Law  for  Isolated  Roughness  -  (Comparison 
of  Present  Investigation  and  Holl  -  1958). 
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region.  However,  significant  differences  may  occur  for  cases  outside 
the  scope  of  the  present  investigation. 

As  mentioned  in  Section  3,4,  Holl  [4]  assumed  that  the  characteristic 
length  ( L^)  is  closely  related  to  the  height  of  the  irregularity  (h) . 

Thus,  it  was  assumed  that  is  equal  to  h.  However,  one  of  the  main 
reasons  for  the  discrepancy  between  the  result  calculated  by  the 
characteristic  velocity  theory  and  the  experimental  data  may  be  due 
to  the  above  assumption.  Therefore,  h  in  equation  (3-10) 
should  be  changed  to  as 


£fp 


L  1 

<-f)m 


£fp 


for 


(5-24) 


and 


°lfp 


(5-25) 


Furthermore,  L  was  derived  semi-empirically  so  that  good  agreement 
c 

between  the  two  results  could  be  obtained. 

Employing  equation  (3-7)  for  the  experimental  value  of 

(Ooe  )  gives 

&fp-exp  ° 


Uc  * 

0?.fp-exp  =  °£fp 


(5-26) 


or 


U  2  0 

(-_£)  , 

-U  ' 


&fp-exp 
°2.fp 


(5-27) 


By  employing  equation  (5-22),  one  can  express  °£fp_eXp 


in  the  form 
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0.53  ,  0.2  1.25 

0 0  -  0.0143  ( ~ )  (— )  (G) 

&fp-exp  6  V 


(5-28) 


or  for  a  zero  pressure  gradient 


2,fp-e.xp 


.  0.53  U  5  0.2 
0.147  (|) 


(5-29) 


where  G  is  equal  to  6.447. 


Substitution  of  equations  (5-23)  and  (5-29)  into  equation  (5-27) 
results  in 


2  ,  0.11  U  6  -0.2 

(S)  =  7 -ooo  (h  (---) 

VU 


(5-30) 


By  knowing  equation  (3-7),  and  by  substitution  of  equation  (5-30) 
into  equation  (5-24)  and  equation  (5-25)  one  getEi 


=  (7.00— -)2  (~) 
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Employing  equation  (3-12)  in  equations  (5-31)  and  (5-32)  results  in 
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(5-34) 


Equations  (5-33)  and  (5-34)  indicate  that  the  dimensionless  character¬ 
istic  length  (L  /6)  for  a  zero  pressure  gradient  is  a  function  of  the 
c 

shape  factor  (H)  and  the  Reynold's  number  bared  on  the  I  oundary  layer 

thickness  in  addition  to  the  dimensionless  relative  height  (h/S) . 

Equation  (5-33)  is  plotted  in  Fig.  5-16  ovei  a  range  of  parameters 

corresponding  to  Fig.  5-8.  It  is  seen  that  in  order  for  L  /6  to  be 

c 

greater  than  one,  h/o  would  have  to  be  significantly  greater  than 
one.  Furthermore,  equations  (5-33)  and  (5-34)  are  plotted  in  Fig.  5~i". 
They  show  that  the  relationship  between  L^/6  and  h/6  is  sensitive  to 
Re^,  especially  for  the  case  where  IW6  is  greater  than  one.  However, 

since  these  two  equations  were  obtained  from  the  results  of  this 

4  5 

investigation  over  the  Re^  range  of  3  x  10  to  2  x  10  and  the  h/6  range 
of  0,06  to  2.50,  care  should  be  taken  in  using  equations  (5-33)  anc 
(5-34)  for  values  of  Re ^  and  h/6  outside  of  the  aforementioned  range 


of  these  parameters. 
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CHAPTER  VI 

SUMMARY,  CONCLUSIONS  and  RECOMMENDATION 


6. 1  Summary 

The  effect  of  pressure  gradient  on  limited  cavitation  was 
experimentally  Investigated  by  using  triangular  protrusions. 

Initially  the  pressure  distribution  was  obtained  in  the  water  tunnel 
test  section  with  the  curved  wall.  Then,  the  velocity  profiles  of 
the  boundary  layer  (6)  on  a  flat  plate  in  the  test  section  were  measured 
by  using  a  total  head  tube  at  seven  stations.  The  nominal  velocity 
was  40  fps. 

The  shape  factor  (G)  and  the  wake  parameter  (IT)  were  calculated 
by  employing  the  pressure  gradient  and  the  boundary  layer  parameters 
with  the  assistance  of  computer  programs  based  on  the  method  of 
E.  Truckenbiodt  [35].  In  so  doing,  the  values  obtained  by  the  method 
of  E.  Truckenbiidt  were  used  for  further  analysis,  after  confirming 
that  there  was  an  insignificant  difference  between  the  boundary  layer 
thickness  obtained  from  the  experiment  and  that  obtained  from  the 
method  of  E.  Truckenbrodt . 

The  limited  cavitation  numbers  were  obtained  visually  employing 
the  desinence  procedure  over  a  range  of  free  stream  velocities  for 
the  various  cases  of  zero,  favorable  and  unfavorable  pressure 
gradients.  The  maximum  range  of  the  velocity  was  30  to  60  fps  and  the 
gas  content  was  maintained  less  than  5  ppm  during  the  experiment. 

Finally,  for  the  case  of  zero  pressure  gradient,  the  effects  on 
limited  cavitation  due  to  the  relative  height  and  the  Reynold’s  number 


based  on  the  boundary  layer  thickness  were  investigated.  The  power- 
law  relationship  between  the  cavitation  number  and  the  aforementioned 
two  parameters  was  obtained  as  in  Ref.  [33].  In  addition,  a 


comparison  of  versus  h/ 6  as  obtained  in  this  investigation  with 

that  of  Ref.  [29],  corrected  to  a  gas  content  of  3.7  ppm  as  described 
in  Appendix  B  was  made  to  obtain  a  more  universal  relationship. 
Furthermore,  the  expression  for  characteristic  length  ( L ) ,  used  in 
the  characteristic  velocity  theory,  was  obtained  semi-empirically 
from  the  above  result. 

For  the  cases  of  favorable  and  unfavorable  pressure  gradients, 
the  shape  factor  (G)  was  found  as  the  best  parameter  to  be  related 
to  pressure  gradient.  The  relationship  for  zero  pressure  gradient 
was  extended  to  include  the  influence  of  shape  factor  (G)  for  a  flow 
with  pressure  gradient. 


6 . 2  Conclusions 

The  following  conclusions  are  based  on  results  from  the  present 
invest igations : 

1.  Hie  limited  cavitation  characteristics  of  isolated  roughness 
for  a  flow  without  pressure  gradient  can  be  represented  in  the  form 


.  a  U  6  b 

U  130 

=  constant  (-j)  (~—  ) 


(6- 


where  a  =*  0.53  and  b  =  0.20  were  obtained  from  the  tests  of  the  two 
dimensional  triangular  protrusions  in  this  investigation. 

2.  The  power  number  'a'  from  the  current  investigation  was  equal 
to  0.53  compared  to  a  value  of  0,40  obtained  in  Ref.  [33],  for  a  flow 


without  pressure  gradient.  The  difference  between  the  two  values 
may  be  due  to  a  difference  in  gas  content. 
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3.  It  was  reaffirmed  that  the  characteristic  velocity  theory 
is  a  rough  estimate  of  0^  for  a  separated  protrusion  element. 

4.  By  employing  a  semi-empirical  expression  for  L^,  used  in  the 

characteristic  velocity  theory,  the  theory  becomes  more  accurate  for 
estimation  for  the  case  of  a  separated  protrusion. 

5.  The  effect  of  pressure  gradient  on  limited  cavitation  for  the 
triangular  protrusion  contradicts  the  characteristic  velocity  theory. 
This  discrepancy  may  be  due  to  disruptions  in  the  local  flow  caused 

by  separation  from  the  protrusion.  Separation  effects  are  not  accounted 
for  in  the  characteristic  velocity  theory. 

6.  The  shape  factor  (C-)  was  toe  most  consistent,  parameter  among 
all  others  considered  to  be  employed  as  a  measure  of  pressure 
gradient  effects  on  the  limited  cavitation  number. 

7.  Equation  (6-1)  was  extended  for  a  flow  with  pressure  gradient, 
to  include  the  influence  of  shape  factor  (G) ,  resulting  in  the  general 
equal  ion 

b  e 


°lfp  "  <4 3  (V>  (G> 


(6-2) 


where  C  =  0.0143,  a  =  0.53,  b  =  0.2  and  e  =  1.25  were  obtained  from 
the  tests  of  the  two  dimensional  triangular  protrusions  in  this 
investigation.  This  correlation  was  shown  to  be  in  good  agreement 
with  experimental  data.  It  is  felt  that  equation  (6-2)  ~an  be  applied 
to  other  types  of  isolated  roughness  by  varying  the  constants  c,  a, 
b  and  e,  which  must  be  obtained  empirically. 
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6 . 3  Recommendat ions 

The  following  is  the  list  of  recommendations: 

1.  Care  should  be  taken  i-  "sing  equation  (6-2)  for  predicting 
the  cavitation  characteristics  of  surface  roughness  at  large  values 
of  the  shape  factor  (G) ,  where  this  relationship  may  not  be  valid. 
Presently,  the  estimations  obtained  with  this  method  are  in  good 
agreement  with  the  data  for  values  of  G  <  7.647. 

2.  Visual  observation  of  desinent  cavitation  should  be  carried 
out  as  far  as  possible  at  constant  gas  content.,  However,  equation 
(5-17)  may  be  useful  for  correcting  for  variations  ir.  gas  content. 

3.  The  cavitation  characteristics  of  additional  roughness 
elements,  including  streamlined  surface  protrusions,  should  be 
invest igated. 

4.  To  apply  the  characteristics  velocity  theory  to  other  types 

of  surface  roughness,  one  must  obtain  for  the  particular  roughness. 

For  a  roughness  with  separated  flow,  roust  be  obtained  experimentally 

whereas  for  the  non-separated  flow  about  a  roughness  one  can  obtain 

0'nc  theoretically. 

Jofp 

5.  As  indicated  under  item  5  in  the  conclusions,  the  discrepancy 
between  experimental  data  and  the  characteristic  velocity  theo  y  may 
be  due  to  separation  effects.  Thus  it  is  recommended  that  experiments 


with  various  pressure  gradients  as  employed  In  this  investigation  be 
conducted  with  a  protrusion  which  produces  a  nonseparated  flow. 


1  ■  -  't*7iw. 
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APPENDIX  A 

DESIGN  OF  EXPERIMENT 


A. I  Curve.- d  Wall 

The  objective  of  mounting  the  curved  wall  in  the  water  tunnel  test 
section  was  to  obtain  various  pressure  gradients  in  the  test  section. 

The  minimum  pressure  coefficient  has  to  be  more  than  -0.4,  derived  by 
the  superposition  equation,  in  order  to  make  a  successful  experiment 
in  the  water  tunnel  because  of  pressure  limitations. 

Originally,  a  NACA  0013-34  airfoil  was  selected  from  reference  [37]. 
Furthermore,  the  minimum  pressure  coefficient  in  this  test  section  was 
derived  by  using  a  computer  program.  Unfortunately,  the  value  was 
found  to  be  below  -0.4  „  so  a  curved  wall  was  modified  to  obtain  a 
proper  value  of  the  minimum  pressure  coefficient.  The  configuration 
and  the  pressure  coefficient  are  shown  in  Fig.  4-5  and  Fig.  5-3, 
respectively. 

A. 2  Triangular  Protrusion 

After  calculating  the  pressure  gradient,  seven  points  were  selected 
for  locating  the  triangular  protrusions.  Three  of  them  were  located 
where  the  pressure  gradient  is  the  maximum,  the  minimum  and  zero.  Four 
other  points  were  selected  such  that  two  of  them  have  the  same  absolute 
pressure  gradient.  The  locations  of  the  seven  stations  are  shown  in 
Fig.  A- 1 . 

The  heights  of  the  protrusions  were  selected  so  that  the  cavitation 
numbers  for  three  cases  within  the  boundary  layer  and  for  one  case 
outside  of  the  boundary  layer  could  be  measured.  The  boundary  layer  was 
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roughly  calculated  by  using  the  equation  for  a  flat  Hate.  The  initial 
schedule  of  the  experiments  is  shown  in  Table  A-l. 
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TABLE  A-l 


INITIAL 

SCHEDULE 

FOR 

CAVITATION 

TEST 

STA. 

I 

II 

III 

IV 

V 

VI 

VII 

<5( inch) 

0.045  0 

1.095  0 

.117 

0.140 

0.167 

0.202 

0.218 

h(inch) 

1/100 

0 

2/100 

0 

3/100 

0 

0 

0 

5/100 

0 

0 

0 

0 

7/100 

0 

0 

0 

0 

0 

0 

0 

10/100 

0 

0 

0 

0 

15/100 

0 

0 

0 

0 

30/100 

0 

0 

0 

0 
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APPENDIX  B 

CORRECTION  OF  CAVITATION  NUMBER  FOR  A  VARIATION  IN  GAS  CONTENT 

The  effect  of  gas  content  for  limited  cavitation  number  has  been 
demonstrated  by  Billet  and  Holl  [38]  as 


°i  ‘  -  cp  ,  + 

min 


kot8' 
1/2  <«£ 


(B-l) 


where  a,  8'  and  k  are  the  dissolved  gas  content  (PPM),  Henry's  law 
constant  (PSI/PPM)  and  an  empirical  adjustment  factor,  respectively. 
Equation  (B-l)  was  initially  proposed  by  Holl  in  Ref.  [29].  For  the 
present  result  and  that  of  Ref.  [29],  equation  (B-l)  becomes 


mxn 


top6’ 

1/2  pii 


or 


a 
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ky 

1/2  Sir 


(B-2) 


(B-3) 


where  the  subscripts  P  and  R  stand  for  present  and  Ref.  [29], 
respectively.  By  using  the  data  of  Ref.  [29],  which  were  carried 
out  under  the  same  condition  except  for  the  gas  content,  k  can  be 
expressed  by  equation  (B-2)  in  the  form 


<oR,  -  V 


(o 


Ra 


°Rb) 


1/2  PU2 

T' 


(B-4) 


where  the  subscripts  a  and  b  are  arbitrary  numbers.  Since  these  two 
experiments  were  carried  out  for  a  water  temperature  of  approximately 
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72°F,  the  Henry's  law  constant  is 

8'  =  1.00  (F SI /PPM)  .  (B-5; 

Substituting  Equation  (B-5)  into  Equation  (B-4)  gives 


k  = 


(0Ra  -  gRb> 

(aF*  -  V 


(B-6) 


The  results  are  tabulated  in  Table  B-l. 

As  indicated  previously,  since  the  present  experiment  was  carried 
out  with  almost  the  same  gas  content,  the  mean  value  in  the  present 
experiment  was  used  as  Thus, 

"  (Vi 

cr  =  — — -  =  3.7  (PPM)  .  (B-7) 

P  n 


The  cavitation  number  of  Ref.  [29],  corrected  for  the  reference  gas 
content  corresponding  to  the  present  investigations,  CJg_corr  is 
expressed  in  the  form  of 


a. 


R-corr 


k(aR  "  V6’ 

1/2  pl£ 


(B-8) 


Employing  equations  (B-5)  and  (B-7)  in  equation  (3-8)  gives 


a 


R-corr 


k(aR  -  3.71) 

1/2  pl£ 


(5-17) 


where  k  is  0. 102 . 


TABLE  B-l 


CONSTANT  (k)  FOR  EFFECT  OF  GAS  CONTENT 


DATA  U  (CtRa  °W 

CO 


SET 

h/6 

(fps) 

(PPM) 

k 

20 

0.048 

1 

0.0634 

30 

10.12 

0.072 

40 

0.088 

50 

0.091 

20 

0.005 

2 

0.2480 

30 

4.49 

0.030 

40 

0.154 

50 

0.117 

20 

0.131 

3 

0.2480 

30 

7.58 

0.166 

40 

0.075 

50 

0.286 

20 

0.028 

4 

0.5030  30 

10.03 

0.067 

40 

0.044 

50 

0.234 

AVERAGE 

0.104 
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